Finite temperature superfluid density in very underdoped cuprates 
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The combination of a large gap, low transition temperature, and quasi two-dimensionality in 
strongly underdoped high temperature superconductors severely constrains the behavior of the ab- 
plane superfluid density p with temperature T. In particular, we argue that the contribution of nodal 
quasiparticles to p{T) is essential to account both for the amplitude of, and the recently observed 
deviations from, the Uemura scaling. A relation between Tc and p(0) which combines the effects of 
quasiparticle excitations at low temperatures and of vortex fluctuations near critical temperature is 
proposed and discussed in light of recent experiments. 



High temperature superconductors appear most strik- 
ingly different from those of standard BCS variety when 
underdoped. While the observed gap in the single- 
particle density of states is large and appears to be only 
increasing with underdoping []|, the a6-plane superfluid 
density at T = 0, /9(0), at the same time is small and 
known to be continuously decreasing 0. The super- 
conducting transition temperature Tc appears to scale 
approximately linearly with yo(0), and this simple rela- 
tion between the two fundamental quantities has been 
ubiquitous in cuprates and widely known as the Uemura 
scaling . Recently, new measurements of the penetra- 
tion depth in the strongly underdoped YBa2Cu307~s 
(YBCO) suggested a significant deviation from the Ue- 
mura relation [3, Q. In this note we point out that 
the combination of three characteristics of underdoped 
cuprates, namely quasi two-dimensionality, a large su- 
perconducting gap, and low T^, irrespective of what the 
microscopic mechanism behind the latter two may be, 
severely constrains the form of piT)- In particular, we 
argue that even the early data of 2] are difficult to ex- 
plain quantitatively by the phase fluctuations alone, and 
that the quasiparticle contribution should make the curve 
Tc{p{0)) convex, as we believe better describes what is 
being experimentally observed Q, 3, 0. 

We begin from the early data on Tc vs. p{0) 0, 
by extracting a quasi-universal number that character- 
izes this relationship 0. The hclicity modulus in a two- 
dimensional (2D) superfluid layer may be expressed in 
units of Kelvin when written as 
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where l^'oP = m*c^/167re^A(T)^ is the superfluid den- 
sity, d a layer thickness, and A(T) the in-plane penetra- 
tion depth. Inserting all the requisite constants one finds 
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In the case of layered quasi two-dimensional materials like 
cuprates there is some ambiguity in defining the length d. 
Assuming it is the minimal thickness along the c-axis so 
that the material still superconducts, d may be estimated 
from the measurements on superconductor-insulator su- 
perlattices or from the crossing-point phenomenon in 



magnetization measurements Using 10 A < d < 12 A, 
with the lower boundary corresponding to the measure- 
ment of Goodrich ct al. ^] on superlattices and the upper 
to the thickness of one c-axis unit cell in YBCO, together 
with Tc = (3.2±0.2)(/im)2K/A2(0) 0, 0, one finds that 
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with the dimensionless amplitude 1.82 < Q < 2.48 in 
YBCO. Similarly, using d = 7.6A and Tc = (3.0 ± 
0.5)(Afm)2K/A2(0) 0, H one finds 1.37 < Q < 1.88 in 
La2-xSrxCuOi (LSCO). Evidently, the amplitude Q is 
not truly universal, but it is a pure number and most 
likely around two. 

First, we argue that assuming that phase fiuctua- 
tions alone reduce the superfiuid density in underdoped 
cuprates 0, ^ leads to difficulties in accounting for 
the observed value of Q. Consider an anisotropic quan- 
tum XY model describing thermal and quantum fluctu- 
ations of the phase of the order parameter, 

S = - J dry^ cos{d{x, z, t) - 9{x + f, z, t)) (4) 
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with X being a discrete two-dimensional (2D) vector la- 
beling sites of a quadratic lattice, 9 ~ drO and r = xi,X2 
a lattice unit vector. The parameter < a < 1 mea- 
sures the c-axis anisotropy. The last term describes the 
quantum fluctuations induced by the interactions and 
to be specific we will assume the Coulomb interaction 
V ^ l/y^\x — + {z ~ z'Y' , although our arguments 
will be rather general. For illustration, we will also con- 
sider the short-range repulsion V = S^^gSz^z'- In general, 
one can write the helicity modulus in the model (4) as 



— = F{ — , — ,a), 
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where F is a dimensionless function of its three dimen- 
sionless arguments and po = pa(J,e,t,a) is the helicity 
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FIG. 1: Helicity modulus vs. temperature in the self- 
consistent Gaussian approximation, for 2D XY model (a = 
e* =0), and with Coulomb interaction (a = 0, = 10, and 
V{q) = The straight dashed line denotes the universal 

BKT limit, p{Tc) = {2/-k)Tc, at which the helicity modulus 
discontinuously vanishes. Note that Tc for fixed po increases 
with interaction. 

modulus at T = 0. (For a = e* = 0, for example, 
Po = ^0 

Let us consider first the 2D limit, a = 0. When e, = 
as well, the last term in (4) completely suppresses quan- 
tum fluctuations, and the Eq. (4) reduces to the stan- 
dard (thermal) 2D XY model. The critical temperature 
for a = e* = is given by Tc = Po/Q2D,xy, with 
Q2D,XY = 1.11 The transition is in the Berezinskii- 
Kosterlitz-Thouless (BKT) universality class, and conse- 
quently at Tc the superfluid densi ty h as a universal dis- 
continuity, with p{T^) = {2/'k)Tc [13 ■ At low tempera- 
tures the reduction of p is linear in T: p{T) ^ po — T/4 
|l2j |. Consider now turning on a weak interaction e* 
while still keeping a = 0. This has two principal ef- 
fects: 1) the reduction of po due to quantum fluctua- 
tions .14]; 2) the introduction of a new (quantum) en- 
ergy scale LUq below which the further reduction of p(T) 
with T is suppressed. In principle, increasing the inter- 
action lifts the energies of the longitudinal (spin-wave) 
phase fluctuations, so Ap{T)/po = {po — p{T))/po should 
be a decreasing function of interactions (see Fig. 1). A 
weak Coulomb interaction e*, nevertheless, is an irrele- 
vant coupling in 2D even at T = "TF| so that the finite 
temperature transition remains in the BKT universality 
class and in particular, p{T~) — {2/-k)Tc still holds for a 
weak e* ^ 0. Defining the dimensionless amplitude for 
e* ^ as Q{a — 0, e*) = Po/Tc, the above discussion 
suggests then that 

2 

- <Q{a^O,e^) <Q2D,XY, (6) 

TT 

which we propose here as a conjecture. 

The lower bound in (6) is dictated by the BKT uni- 
versality class of the finite temperature transition in 2D. 
Our conjecture of the existence of the upper bound on Q 




I , I , I I I I L 

0.2 04 06 08 

2 



FIG. 2: The amplitude Q{a = 0,e.) = po/Tc in the self- 
consistent Gaussian approximation for the helicity modulus 
in 2D, Eq. (7), with Tc defined by p(Tc) = (2/7r)T, for a 
long-range and a short-range interactions. The lines are the 
weak-interaction results in Eq. (8), and J = 1. 

in two dimensions is supported by a self-consistent cal- 
culation in which one approximates the first term in Eq. 
(4) (with a = 0) by an optimally chosen Gaussian term 
[Tl| . The straightforward calculation then gives 

where F^{q) — V{q)^fSin'^{^). Identifying the ap- 
proximate transition temperature as the point when clas- 
sical vortices unbind and the self-consistent Gaussian ap- 
proximation breaks down, p(Tc) = (2/t:)Tc, at weak in- 
teraction Eq. (7) yields 

Q(a = 0,e,) = — [1--^ / ^F{q) (8) 
71" 4J2 J [2ny 

Numerical results for Q{a = 0, e*) in our self-consistent 
Gaussian approximation are shown in Fig. 2. In both the 
Coulomb and the short-range case we find Q to be a de- 
creasing function of interaction. Although the above cal- 
culation in principle includes only the longitudinal (spin- 
wave) phase fluctuations, it should be a reasonable indi- 
cator of the interaction dependence of the amplitude Q. 
Recall that the main contribution to A/9(T) at temper- 
atures T < 0.8Tc in the classical 2D XY model actually 
comes from the spin- waves [T^ . and it is the increase of 
their energies with interaction that should ultimately be 
responsible for the proposed behavior of Q(a = 0, e,). 

At low temperatures T <^ ujq, from Eq. (7) one also 
finds that for Coulomb interactions in 2D 

^ = l-const.( — )^ (9) 

Po 
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where ujq = et^^PQ^^ is the characteristic temperature 
scale for quantum fluctuations. For a short-range inter- 
action one similarly finds 1 — {p{T)/ po) ^ {T/ujq)^ in 2D, 

2/3 2/3 

with ujq = ej Pq . Both expressions show the antici- 
pated reduction of Ap{T)/po with interactions for fixed 
Pq. Also, note that when pQ ^ 0, uiq ^ po in general in 
2D. This will have important consequences for the form 
of p{T), as we will discus shortly. 

The lesson from the above considerations is that in- 
teractions are expected to reduce the amplitude Q from 
its value in the 2D XY model. Introducing the third 
dimension by turning on a > in (4) does the same. 
This is evident in the Monte Carlo results ^3 "^^11 
as from the 3D value Q3d,xy = 0.45 0, for example. 
The reason behind this trend is quite physical: coupling 
the layers only strengthens the phase coherence so that 
a higher temperature is required to destroy it. Both the 
interactions and the coupling between the layers should 
decrease Q and we expect that 

Q(a,e*) < Q2Axy = 1.11 (10) 

for the family of quantum XY models defined by Eq. (4). 
Introducing longer-ranged couplings than just between 
the nearest neighbors as long as there is no frus- 
tration in the problem 20], would also increase the ten- 
dency to order and consequently should reduce Q further. 
The maximum value of Q2D.XY then corresponds to the 
weakest (unfrustrated) superfluid order on the quadratic 
lattice [2IJ. 

The conjectured maximum value of the ratio Q in the 
family of quantum XY models (4) appears to be signifi- 
cantly below the one characterizing the standard exper- 
imental Uemura plot. To bring the two closer one needs 
to assume a considerably smaller 'thickness' d: choos- 
ing the average spacing between the adjacent layers in 
YBCO, d = 6A, for example, would lead to Q > 1.09 
from experiment, and possibly falling just slightly below 
our upper bound. A more natural interpretation of the 
above discrepancy is that an additional mechanism for 
the reduction of superfiuid density is required. Of course, 
one is readily available in cuprates: nodal quasiparticles. 

As well known [2^ |2^ . with the phase fluctuations 
neglected, at low temperatures depletion of the helicity 
modulus in a 2D d-wave superconductor is entirely due 
to low-energy quasiparticles near the four nodes of the 
order parameter, 

p = p{0)-yT + O{^). (11) 

The number y = {2ln{2)/Tr)z'^{vp /va), A is the max- 
imum of the superconducting gap, vp and va sue the 
two velocities characterizing the low-energy quasiparti- 
cle dispersion, and z is the charge renormalization fac- 
tor. Assuming A ;2> Tc in underdoped cuprates makes 
the correction to linear dependence in Eq. (7) negligible 
for T < Tc and we drop it hereafter. 

Very underdoped cuprates are characterized by a large 
gap A ^ Tc, small transition temperature Tc, and strong 
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FIG. 3: Schematic behavior of the superfluid density vs. tem- 
perature in a low- Tc, large gap, quasi-2D d-wave superconduc- 
tor. Doping is assumed to be decreasing from top to bottom. 
Note that both Tc and the slope should decrease with un- 
derdoping (see the text). Dashed line again represents the 
universal BKT limit. 

c-axis anisotropy. We may then include both the quasi- 
particle and phase fluctuation contributions to the re- 
duction of the superfluid density by reasoning as follows: 
1) Assume that the T = superconducting transition in 
underdoped regime is continuous. General scaling argu- 
ments imply then that near the quantum critical point 
p{0) - x(D+z-2)/z^ ^^.^^ for £) = 2 (2D) reduces to 
p(0) ^ Tc independently of the nature of the putative 
quantum critical point. It then follows that as Tc —>■ 0, 
the quantum temperature scale in Eq. (9) ujq ^ Tc 
and the longitudinal phase fluctuations are strongly sup- 
pressed for all T < Tc- 2) With phase fluctuations being 
essentially negligible, quasiparticles alone reduce p{T) ac- 
cording to Eq. (11), with only the linear term being im- 
portant for T < Tc <C A. 3) Finally, the superfiuid den- 
sity in 2D still can not fall below the universal BKT value, 
since for p{T) < {2/-k)T the superfiuid order would be- 
come unstable to unbinding of classical vortices. In sum, 
the transition temperature in a quasi-2D d-wave super- 
conductor with A Tc is approximately determined by 

IT 

p{Q)-yTc^—, (12) 

TT 

where p{Q) is the physical T — helicity modulus, already 
renormalized by quantum fluctuations. p{T) at various 
dopings acquires then an approximate simple form de- 
picted schematically in Fig. 3. Virtual (bound) vortex 
fluctuations should additionally reduce p{T) near Tc, and 
thus tend to slightly lower the transition temperature. 
With < a <C 1, on the other hand, some rounding of 
the universal BKT jump close Tc is expected, as well as 
a slight increase of Tc- These fine effects are being omit- 
ted from Fig. 3 for simplicity but should be expected in 
experiments. 

The addition of the quasiparticle contribution y in Eq. 
(12) makes the agreement with old Uemura data imme- 
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FIG. 4: Data reported in Q- The lines are fits to the data as 
discussed in the text. The convexity of the curves is due to 
the doping dependence of the quantity z^{vf/va)- 

diately better: taking y ~ 1 near optimal doping p3 | 
brings the sum in Q ~ i/ + (2/7r) closer to two. More 
importantly, Eq. (12) implies a systematic deviation of 
Tc(p(0)) from the simple linear relation due to the antic- 
ipated decrease of the slope y with underdoping. Even if 
one assumes the charge renormalization factor z to ap- 
proximately stay constant with doping, the decrease of 
VA Q together with the approximate constancy of vp 
[25l makes the slope y inevitably an increasing function 
of doping, and therefore of p(0). The slope of the Tc(p(0)) 
curve should therefore be a decreasing function of p(0), i. 
e. the curve should be convex upward. In fact we expect 
that the dependence of y on doping should be stronger 
than what would follow from the observed increase of va 
alone. This would be in agreement with the general the- 
ory of the fluctuating d-wave superconductors [26j and 
the gauge theory of t-J model [23, which both predict 
the charge renormalization factor z oc Tc at small dop- 
ings. 

Figure 4 shows the data taken from Pereg-Barnea et 
al. Q where circles (squares) are from a(6)-axis mea- 
surements. To avoid contributions arising from chains in 



the 6-direction, we consider only the a-direction, though 
both data sets clearly show a departure from linear be- 
havior. In fitting the data we used a straight line fit to 
the doping dependence of the superconducting gap Aq 
from Sutherland et al. (their Figure 6) [l| and converted 

to VA using Aq — hkpVA/'^ with kp ~ O.TA . We have 
taken vp = 2.5 x 10'' cm/s. The dotted line represents 
then the Eq. 12 with d = lOA and z = 1 and no fitting 
parameters. A somewhat better agreement is obtained 
by assuming z = Tc/C with C = 79K, which we show as 
the full line, as an illustration. The reader is cautioned, 
however, that such a strong dependence of the slope y on 
Tc has not yet been directly observed. The dashed line 
is the linear Uemura relation with Q — 2, demonstrating 
that the data departs significantly from this behavior. 

The curvature in the Tc(p(0)) plot was also noted and 
discussed in a recent paper by Tallon et al. Q where it 
was attributed to the opening of an additional gap in the 
density of states competing with superconductivity. Our 
explanation is fundamentally different in nature and in 
particular, the curvature is argued to be a generic fea- 
ture of two-dimensional d-wave superconductor with the 
increasing gap and the decreasing Tc. The microscopic 
mechanism behind the latter two properties of under- 
doped cuprates is left unaddressed here and will be a 
subject of a separate publication [2^ . 

In sum, we discussed the reduction of superfluid den- 
sity in a quasi-2D d-wave superconductors with low Tc 
and large superconducting gap. We argued that the 
Coulomb interaction under these conditions essentially 
eliminates the longitudinal phase fluctuations but that 
the transition at Tc should still be governed by the uni- 
versal vortex physics and fall into the BKT universality 
class. A relation between p(0) and Tc which takes into 
account both the quasiparticle contribution and the uni- 
versal vortex fluctuations near Tc is proposed to explain 
the deviation from the linear Uemura scaling observed in 
recent experiments on underdoped YBCO. 
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